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We analyze the physics of bright solitons in two-dimensional dipolar Bose-Einstein condensates. 
These solitons are not possible in short-range interacting gases. In particular, we discuss the neces- 
sary conditions for the existence of stable 2D bright solitary waves. The stability of the solutions is 
studied by means of the analysis of lowest-lying excitations. Additionally, we study the scattering of 
solitary waves, which contrary to the contact-interacting case, is inelastic, and could lead to fusion 
of solitary waves. Finally, the experimental possibilities for observability are discussed. 



During the last years, the physics of ultracold atomic 
and molecular gases has attracted a large interest. Al- 
though these gases are very dilute, their properties are 
crucially determined by the interparticle interactions . 
Up to now, only short-range van der Waals interactions 
have played a significant role in typical experiments. 
However, very recent developments are paving the way 
towards a new fascinating research area, namely that of 
degenerate dipolar gases. A major breakthrough has 
been very recently performed at Stuttgart University, 
where a Bose-Einstein condensate (BEC) of Chromium 
atoms has been realized .2]. Chromium atoms are par- 
ticularly interesting, since they present a large magnetic 
dipole moment, fi — 6/ib- Hence, Chromium BEC con- 
stitutes the first realization of a degenerate dipolar gas. 
On the other hand, recent developments on cooling and 
trapping of molecules Q, on photoassociation experi- 
ments!^, and on Feshbach resonances of binary mix- 
tures |3, open exciting perspectives towards the realiza- 
tion of a degenerate gas of polar molecules. Particularly 
interesting in this sense is the recent creation of RbCs 
molecules in their absolute vibronic ground state X 1 S + , 
with a calculated electric dipole moment of 1.3 D jfj. 

In addition to the usual short-range forces, dipolar par- 
ticles oriented by some external electric or magnetic field 
interact via dipole-dipole interaction, which is long-range 
and anisotropic, being partially attractive. New exciting 
physics is therefore expected in these systems. In par- 
ticular, recent theoretical analyses have shown that the 
stability and excitations of dipolargases are crucially de- 
termined by the trap geometry pjlalSi- Ultracold dipolar 
particles are also attractive in the context of strongly- 
correlated atomic gases as ph ysic al implemen- 
tation of quantum computation ideas [12J, and for the 
study of ultracold chemistry 0] . 

The nonlinearity of the BEC physics is one of the ma- 
jor consequences of the interparticle interactions. In this 
sense, resemblances between BEC physics and nonlin- 
ear physics (in particular nonlinear optics) have been an- 
alyzed in detail. Several remarkable experiments have 
been reported in this context, including four- wave mix- 
ing 0|> BEC collapse Fiji , and the creation of dark, 
bright and gap solitons 16] . The physics of BEC soli- 



tons has indeed arouse a large interest. For short-range 
interactions and at sufficiently low temperature, the BEC 
physics is provided by a nonlinear Schrodinger equation 
(NLSE) with cubic nonlinearity (Gross-Pitaevskii equa- 
tion) [l|. In ID this equation admits solitonic solutions 
for bright (dark) solitons [17| for attractive (repulsive) 
interatomic interactions, the equivalent of self-focusing 
(self-defocusing) nonlinearity in Kerr media. 

The equation determining the physics of a dipolar BEC 
is, as we discuss below, a NLSE with nonlocal nonlinear- 
ity, induced by the dipole-dipole interaction. Some in- 
teresting effects related with nonlocal nonlinearity have 
been studied in different contexts. In particular, it has 
been suggested that the nonlocality may play an impor- 
tant role in the physics of solitons and modulation insta- 
bility . Recent experiments in nematic liquid crystals 
have confirmed some of these predictions 19]. Particu- 
larly interesting is the possibility of stabilization of local- 
ized waves in cubic nonlinear materials with a symmetric 
nonlocal nonlinear response [2oj] . 

In this letter we are interested in the physics of soli- 
tons in two-dimensional dipolar condensates. Whereas 
for short-range interacting BEC a 2D stable solitary wave 
is not possible, we show that the dipole-dipole interac- 
tions may stabilize a 2D solitary wave. We discuss the 
appropriate conditions at which this is possible, which 
involve the tun ing of the dipole-dipole interactions as dis- 
cussed in Ref. |2lJ . We study the stability of the solitary 
waves, discussing its lowest-lying excitations. The scat- 
tering of 2D solitary waves is also analyzed. We show that 
this scattering is inelastic, and it can lead under appro- 
priate conditions to the fusion of solitary waves. Finally, 
we conclude with a discussion about the experimental 
feasibility and generation of the solitary waves. 

In the following we consider a BEC of N particles with 
electric dipole d (the results are equally valid for mag- 
netic dipoles). The dipoles are oriented by a sufficiently 
large external field, and hence interact via a dipole-dipole 
potential: 

V d {r) = 5d (l-3cos 2 0)/r 3 , (1) 

where gd = aNd 2 /Aireo, with eo the permittivity of 
vacuum, 8 is the angle formed by the vector joining 
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the interacting particles and the dipole direction, and 
a is a parameter tunable as follows. As in Ref. [2l| 
we consider that the electric field which orients the 
dipoles is a combination of a static field along the z- 
direction and a fast rotating field in the radial plane 
E{t) — E{cos4>z + sm0[cos(fM)i + sin(flt)y}}. The ro- 
tating frequency Q is sufficiently slow to guarantee that 
the dipole moments follow adiabatically the field E(t), 
but sufficiently fast such that the atoms do not signifi- 
cantly move during the time j2l|. In this case, the 
interaction can be averaged over a rotation period, re- 
sulting in a tunable dipole-dipole potential, with a = 
(3 cos 2 4> — l)/2, which may range from 1 to —1/2. The 
tunability of the dipole-dipole interaction will become 
crucial in our discussion of 2D solitary waves. 

The physics of the dipolar BEC at sufficiently low tem- 
peratures is provided by a NLSE with nonlocal nonlin- 
earity 0: 
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where J \ijj(r,t)\ 2 dr = 1, and g — Anh 2 aN/m is the cou- 
pling constant which characterizes the contact interac- 
tion, with a the s-wave scattering length. In the fol- 
lowing we consider a > 0, i.e. repulsive short-range in- 
teractions. We assume an external trapping potential 
U(f) = muj 2 z 2 /2, with no trapping in the xy-planc. 

In the following we are interested in the possibility to 
achieve a 2D localized wave in dipolar BECs. In order to 
have a good insight on this issue, we introduce a Gaussian 
Ansatz for the wavefunction: 



1 



tVhI /2 l p lY- 



■ exp 



V 



2l 2 L 2 



(3) 

where l z = y^H/muj z , and L p and L z are dimensionless 
variational parameters related with the widths in the xy- 
plane and the z-direction, respectively. Using this Ansatz 
the energy of the system reads: 
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where g = 2gjhio z l\ — 8ttNo,/I z and = 2g c i/fuJ z l^ : , and 
f(n) = (k 2 - I)" 1 [2k 2 + 1 - 3k 2 H(k)] , with H(n) = 
arctan(V k 2 — 1) / \/k 2 — 1) . 

The minimization of E leads to the equations: 
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where F(k) = (1 - n 2 )- 2 [-4k 4 - 7k 2 + 2 + 9k 4 F(k)] , 
G(k) = (1 - k 2 )" 2 [-2k 4 - +10k 2 + 1 - 9k 2 H (k)] , and 
(3 = cjd/g- Eqs. J5J and JHJ admit a solution, and hence 
a localized wave, only under certain conditions. A sim- 
plified picture may be achieved by considering the fully 
two-dimensional situation in which the confinement in z 
is strong enough to guarantee L z = 1. In that case, both 
kinetic and interaction energy scale as l/L 2 . In absence 
of dipole-dipole interactions (gd — 0), and irrespective 
of the value of L p , E(L p ) is always either growing with 
L p (collapse instability) or decreasing with L p (expan- 
sion instability). This reflects the well-known fact that 
2D solitons are not stable in NLSE with contact inter- 
actions. In the case of a dipolar BEC, the situation is 
remarkably different, since the function / depends explic- 
itly on L p . This allows for the appearance of a minimum 
in E(Lp) (inset Fig.^J, which from the asymptotic values 
of / (/(0) = —1 and /(k — » oo) = 2) should occur if: 
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A simple inspection shows that this condition can be ful- 
filled only if <jd < 0, i.e. only if the dipole is tuned as 
previously discussed with tf> > 54.7° (this is true also for 
L z 7^ 1). In that case, the tuning of the dipole-dipole 
interaction may allow for the observation of a stable 2D 
solitary wave, characterized by an internal energy Eg < 
(inset Fig.^). Note that if Na/l z 3> 1, then we arrive to 
the condition \j3\ > 3/8tt ~ 0.12. We have compared this 
value with the result obtained from the direct resolution 
of Eq. (|10fl (see below) for the case of large g, obtain- 
ing stable 2D solitary waves for \(3\ > 0.12, in excellent 
agreement with the Gaussian Ansatz. 

In order to gain more understanding on the stability 
of the 2D solitary waves, we have analyzed the lowest ly- 
ing modes of these structures, namely the breathing and 
quadrupolar modes. To this aim, we employ a Gaussian 
Ansatz of the form El El: 



*(f;<) = * ^, ^] e^ x2+ ^y 2+ ^y 2 (8) 
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where {&»(£), Pi{t)}, i — x,y,z, are time-dependent pa- 
rameters, and insert this Ansatz into the corresponding 
Lagrangian density 

C = —^y*-i,y*) + —\vy\2 + l muJ l z 2\y(r 7 t)\ 2 + 
2 2m 2 
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After integrating L = J drC, we obtain the correspond- 
ing Euler-Lagrange equations for {bi(i), fli{t)}. Lineariz- 
ing these equations around the stationary solution ob- 
tained from Eqs. © and ©, we obtain the expressions 
for the frequencies of the lowest-lying modes [23j . A typ- 
ical behavior of the frequencies of the lowest modes is 
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depicted in Fig. The lowest-lying mode has for any 
value of P a breathing character. For sufficiently small 
values of the breathing mode tends to zero, and even- 
tually the system becomes unstable against expansion. 
This corresponds to the disappearance of the minimum 
in the inset of Fig.^ In this regime, the 2D picture pro- 
vides a good description of the physics of the problem, as 
shown in Fig. ^ For sufficiently large values of \j3\, the 
3D character of the system becomes crucial, leading to 
a different sort of instability, in this case against 3D col- 
lapse. This is reflected in the decrease of the frequency 
of the breathing mode. 
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Figure 1: Breathing (bold dashed) and m = ±2 quadmpole 
(bold solid) mode for g = 20, with ko = Lp/L z . Results from 
a purely 2D calculation are shown in thin lines. Inset: E(k) 
for g = 500, and /3 = -0.10 (dashed) and (3 = -0.20 (solid). 

As shown in Fig. ^ a 2D calculation offers a good 
description of the problem for sufficiently small val- 
ues of |/3 1 . In that case the wavefunction factorizes 
as = ip(p)ipo(z), where ipo is the the ground- 

state of the harmonic oscillator in the z-direction. Em- 
ploying this factorization, the convolution theorem, and 
the Fourier transform of the dipole potential, Vd(k) = 
(47r/3)(3fc 2 /fc 2 - 1), we arrive at the 2D NLSE: 
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where h is the Fourier transform of n(p) = |V>(/?)| 2 , and 
h<2D(k) = — 2 — 3y/irke k erfc(fc), with erfc(x) the comple- 
mentary error function. Below, we employ Eq. (|10ll to 
analyze numerically the dynamics of 2D solitary waves. 

Up to now, we have analyzed a single localized wave, 
showing that a stable solitary wave may exist under ap- 
propriate conditions. In order to deepen our understand- 
ing of these 2D solutions, and their comparison with the 
solitonic solutions of the ID NLSE, we have analyzed 
the scattering of two of these localized waves for differ- 
ent values of their initial center-of-mass kinetic energy, 



Ekin- Direct numerical simulations of the 2D nonlocal 
NLSE show that the scattering of dipolar 2D solitary 
waves is inelastic. In particular, as shown in Fig. [21 for 
sufficiently slow localized waves (for the case considered 
in Fig. |21 Ekin < 2.9|i?s|) two solitary waves merge when 
colliding. As observed in Fig. [21 the solitary waves, when 
approaching, transfer their center-of-mass kinetic energy 
into internal energy, transforming into a single localized 
structure. This structure, although localized, is in an ex- 
cited state, and clear oscillations may be observed. For 
larger initial kinetic energies, the waves move apart from 
each other after the collision, but the transfer of kinetic 
energy into internal energy is enough to unbind the soli- 
tary wave, and the solitary waves are destroyed. This in- 
elastic character of the scattering of dipolar 2D solitary 
waves, clearly differs these solutions from the solitonic 
solutions of the ID NLSE ^t|> an d will be the subject of 
further investigation. 
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Figure 2: Density plot of the fusion of two dipolar 2D solitary 
waves for g = 20, /3 = —0.5, and k = 0.01. From top to 
bottom co z t/2 = 0, 1000, 2000, 3000, 4000, 5000. 

In the final part of this Letter, we would like to dis- 
cuss some issues concerning the experimental realization 
of stable 2D solitary waves in ultracold dipolar gases. 
As we have previously mentioned, the experimental gen- 
eration of these structures demands the tuning of the 
dipole-dipole interaction, and in particular the inversion 
of its sign. In addition to this condition, the dipole-dipole 
interaction must be sufficiently large, \gd\/g > 0.12. Un- 
fortunately, for the case of Chromium, \gd\/g — 0.03 for 
a tuning angle 4> = 7r/2, and hence the tuning of the 
dipole-dipole interaction must be combined with a reduc- 
tion of the contact interactions via Feshbach resonances. 
This combination results problematic due to technical 
limitations Therefore, the best candidate for the 
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generation of 2D solitary waves is a condensate of polar 
molecules. As commented in the introduction of this Let- 
ter current developments open very optimistic perspec- 
tives for the achievement of a BEC of polar molecules in 
the very next future. Since appropriate molecules in the 
lowest vibronic state can present very large electric dipole 
moments [(J, and the dipole-dipole interaction may be 
tuned in these gases by means of rotating electric fields, 
a broad regime of values of f3 may be available, and hence 
a stable 2D solitary wave should be observable. 




x/l z 

Figure 3: Density profile for g — 200, and f3 — —0.10 (solid) 
and —0.20 (dashed). For both cases a xy trap with oscillator 
length l± = 401 z has been employed. In the inset, width of 
the BEC as a function of |/3| for the same parameters. 

In the previous discussion we did not consider any trap- 
ping on the xy plane. In Fig. |3| we consider the case of 
a loose xy-trapping. A large shrinking of the BEC den- 
sity profile evidences the appearance of the solitary-wave 
solution. We would like to stress, that this large macro- 
scopic effect may be scanned using the tuning techniques 
in polar molecules, it does not involve any collapse of the 
wavefunction (which is prevented by the dipole-dipole 
interactions), and it occurs for values of |/3| <C 1 for 
which the dipole-dipole interaction could be considered a 
perturbation to the short-range one. Finally, we should 
point out that analysis of the unperturbed propagation of 
the solitary wave may be performed by inducing bounc- 
ings within loose xy-traps, or by adiabatically removing 
the xy confinement. Collisions of solitary waves could 
be performed by generating the waves in two separated 
dipolar traps, and subsequently releasing them. 

In summary, we have shown that contrary to the case 
of short-range interacting BECs, stable 2D solitary waves 
can be generated in dipolar BEC (especially in molecular 
BEC) by means of tuning techniques, if the dipole-dipole 
interaction is sufficiently large. These solitary waves scat- 
ter inelastically, and may undergo fusion for sufficiently 
low scattering energies. For trapped gases, a significant 
shrinking of the BEC is predicted even for ly^l <C g. 
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